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基础 研究 


SR 测度 的 动态 心肌 PET 图 像 因子 分 析 
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Se uM m TU T 
组 织 的 时 间 活 度 曲线 。 方 法 通过 最 小 化 动态 图 像 与 因子 模型 间 的 c 散 度 将 动态 PET 图 像 做 初步 的 分 解 , 得 到 初始 因子 与 因子 
图 像 ,然后 联合 PET 像 素 动力 学 聚 类 的 先 验 信息 解决 因子 分 析 模 型 中 解 的 非 唯一 性 问题 ,将 初始 因子 与 因子 图 像 通过 空间 变换 
生成 具有 生理 意义 的 组 织 活 度 曲线 及 组 织 空间 分 布 。 结 果 与 传统 的 最 小 二 乘法 测度 和 最 小 化 因子 图 像 间 重 苹 程度 约束 模型 
相 比 ,本 模型 对 噪声 的 敏感 性 较 低 ,提取 出 的 结果 的 精确 性 较 高 。 结 论 通过 选取 最 优 的 a 值 作为 因子 分 析 模 型 的 测度 ,并 引入 
PET 图 像 像 素 的 动力 学 聚 类 信息 ,能 精确 地 获得 血 输入 函数 及 局 部 组 织 的 时 间 活 度 曲 线 , 在 视觉 评价 及 量化 评价 均 具 有 优 
质 表 现 。 
关键 词 :因子 分 析 ;a 散 度 ; 正 电 子 发 射 计算 机 断层 显 像 ; 动 力学 此 类 


Kinetic cluster and a-divergence-based dynamic myocardial factorial analysis of positron- 
emission computed tomography images 
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Guangzhou 510515, China 


Abstract: Objective We purpose a novel factor analysis method based on kinetic cluster and a-divergence measure for 
extracting the blood input function and the time-activity curve of the regional tissue from dynamic myocardial positron 
emission computed tomography(PET) images. Methods Dynamic PET images were decomposed into initial factors and factor 
images by minimizing the a-divergence between the factor model and actual image data. The kinetic clustering as a priori 
constraint was then incorporated into the model to solve the nonuniqueness problem, and the tissue time-activity curves and 
the tissue space distributions with physiological significance were generated. Results The model was applied to the “RbPET 
myocardial perfusion simulation data and compared with the traditional model-based least squares measure and the minimal 
spatial overlap constraint. The experimental results showed that the proposed model performed better than the traditional 
model in terms of both accuracy and sensitivity. Conclusion This method can select the optimal measure by a value, and 
incorporate the prior information of the kinetic clustering of PET image pixels to obtain the accurate time-activity curves of the 
tissue, which has shown good performance in visual evaluation and quantitative evaluation. 
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无 创 测定 局 部 心肌 血 流 量 (RMBF) 一 直 是 临床 心 析 , 目 前 广泛 采用 的 估计 血 输 入 函数 方法 为 基于 感 兴 
脏 学 所 追求 的 日 标 之 一 "”。 通 过 从 动态 图 像 中 估计 血 、 ”区 域 法 ,该 方法 简单 易 行 ,但 极度 依赖 于 医生 手工 勾画 
输入 函数 与 输出 函数 后 ,并 应 用 动力 学 模型 ,可 以 得 到 ” ”的 感 兴趣 区 域 ,而 感 兴趣 区 域 的 准确 性 又 受 医 生 个 人 经 
绝对 的 定量 参数 (如 物质 转运 速率 .代谢 速度 和 受 体 验 及 部 分 容积 效应 影响 。 因 子 分 析 作 为 从 变量 群 中 提 
结合 率 等 功能 参数 等 ), 从 而 实现 对 心肌 血 流 的 定量 分 。 取 公 共 因 子 的 统计 技术 ,已 经 被 用 来 从 动态 图 像 中 无 
析 “”。 在 动力 学 模型 的 应 用 中 , 血 输 入 函数 的 估计 是 非 ” 创 地 提取 血 输 入 函数 和 组 织 的 时 间 活 度 曲 线 ” ,如 基 
党 关键 的 ,结果 的 准确 与 否 直接 影响 到 后 期 的 定量 分 ”于 主 成 分 分 析 王 2 .独立 成 分 分 析 “ .最 小 二 乘法 “号 、 
TTE 最 大 似 然 法 “的 因子 分 析 模型 相继 被 提出 。 然 而 ,用 因 
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斯 噪声 模型 的 最 小 二 乘 测度 与 基于 泊 松 噪声 模型 的 最 
大 似 然 测度 “*””。 而 临床 获取 动态 PET 图像 噪声 很 难 
直接 用 高 斯 或 泊 松 模型 进行 刻画 ,因此 ,本 文 提 出 基于 a 
散 度 的 因子 分 模型 ,a 散 度 由 S.Amari 提 出 ,是 诸多 信息 
散 度 族 中 的 一 种 ,通过 选取 不 同 的 wx 值 ,可 以 用 于 任意 噪 
声 模型 的 数据 测度 2 。 

针对 因子 分 析 模型 中 解 的 非 唯一 性 ,传统 上 利用 最 
小 化 因子 图 像 间 的 重 秋 程度 这 一 约束 来 解决 该 问题 ,但 
由 于 受到 PET 图 像 分 状 率 和 部 分 容积 效应 的 影响 ,心肌 
信号 经 常会 摊 杂 着 10%~15% 的 血液 信号 ,基于 最 小 化 
子 图 像 之 间 的 重合 程度 的 约束 未 能 考虑 到 这 一 事 
实 。 在 动态 PET 人 研究 中 ,每 一 个 像素 的 时 间 行 为 可 以 通 
过 曲线 时 间 活 力 (TAC) 来 描述 ,我 们 将 具有 相同 时 间 行 
为 的 像素 进行 聚 类 ” ,把 这 一 聚 类 结果 作为 先 验 约 束 用 
来 解决 因子 分 析 模 型 中 的 非 唯一 性 ,得 到 具有 生理 意义 
的 组 织 活 度 曲 线 。 因 子 分 析 模 型 应 用 于 基于 XCAT 
(Extended Cardiac- Torso phantom) Ji ^* [Jj Cf "Rb 
PET 心 肌 灌注 仿真 数据 中 ,实验 均 表 明 我 们 的 方法 在 视 
党 评价 及 量化 评价 等 方面 优 于 传统 的 因子 分 析 模 型 。 


1 方法 
1.1 PET 图 像 的 获取 

对 于 动态 序列 PET 图 像 忆 像素 点 /在 第 六 帧 的 示 踪 
剂 活 度 zw 可 通过 以 下 积分 式 计 算得 到 交 ， 


Dynamic sequence 


1 对 动态 序列 做 因子 分 析 示 意图 


Factor images 


Xin = fe. t)e"dt (1) 
HP eG, ORIRE Eit WISSEN 
fL tn 2] HERZ E m ERARA E. AA 
Xu Xa t Xin, 
pe cx 2) 
x; X, Hind Xn, 
其 中 ren" 7 ,7 的 一 列表 示 一 帧 图 像 ,7 的 一 行 表 
示 单 个 像素 点 的 时 间 活 度 曲线 。 


动态 PET 序 列 图 像 可 分 解 为 有 限 的 时 间 序 列 向 量 
FE 和 与 其 相对 应 体 素 的 空间 分 布 C, 这 一 分 解 模型 可 用 
下 式 来 表示 

I=CF+e (3) 

其 中 动态 PET 图 像 1 的 大 小 为 Nx7, 其 中 入 表示 每 
一 帧 PET 图 像 的 像素 点 数目 ,7 表示 动态 网 像 的 帧 数 。 
因子 矩阵 的 大 小 为 玉 x7, 因 子 图 像 C 的 大 小 为 WxK， 
其 中 天 为 所 需 提取 的 因子 数目 。s 表 示 剩 余 信 号 , 即 图 像 
的 噪声 和 被 模型 所 忽略 的 小 部 分 信号 。 在 这 里 ,为 保证 模 
型 的 实际 应 用 价值 ,因子 数目 K 要 满足 K<<N,K<<7。 

图 1 表示 对 动态 序列 PET 图 像 做 因子 分 析 的 示意 
图 ,在 理想 状况 下 ,因子 矩阵 正 对 应 于 动态 图 像 中 各 个 
生理 组 织 的 时 间 活 度 曲线 ,因子 图 像 C 对 应 于 各 个 生理 
组 织 的 空间 分 布 。 


-一 Blood 
- - Myocardium 


Activity (n Ci/cc) 


Time (s) 


Factors 


Fig.1 A schematic illustration of factor analysis for dynamic sequences. 


12 提出 的 因子 分 析 模 型 

子 分 析 模 型 的 求解 主要 由 两 个 步骤 来 完成 :(]) 
通过 最 小 化 模型 与 真实 图 像 之 间 的 测度 来 分 解 图 像 , 得 到 
初始 的 因子 与 因子 图 像 。(2) 解 决 模型 中 的 非 唯一 问题 。 


我 们 通过 a 散 度 来 测定 图 像 7 与 模型 CF 之 间 的 偏 
差 ,a 散 度 由 S.Amari 提 出 ,是 诸多 信息 散 度 族 中 的 一 
种 ,用 于 测量 两 个 数据 分 布 p(x) 和 g(x) 之 间 的 偏差 ,其 
定义 为 : 
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DAP) = g gP 0 - 9, p/4/ ") ae R (4) 


结合 信息 散 度 的 特点 ,本 文 提出 最 小 化 基于 a 散 度 
测度 的 目标 函数 来 估计 因子 与 因子 图 像 C: 
al, +(1 -OCP), 
iuret aE en IL EA 90 
基于 a 散 度 的 动态 PET 图 像 的 因子 分 析 模 型 可 表 
示 为 以 下 上 是 优化 问题 : 


(ar mare mum fa Cop) (6) 

我 们 采用 一 种 迭代 的 方式 去 解决 上 式 , 具 体 算法 
如 下 : 

i 初始 化 C,F: 设 置 C 中 元 素 为 0-1 之 间 的 随机 数 ， 
再 根据 C 估 计 初 始 F; 

iif: 


(a)C-step: 根 据 当前 的 动态 图 像 T 和 Fr 估计 CC; 
(b) 将 C 中 的 小 于 零 的 元 素 置 为 0:; 

(c)F-step: 根 据 当前 的 动态 图 像 T 和 CC 估计; 
(d) 将 正中 的 小 于 零 的 元 素 置 为 0:; 

(e) 检 查 是 否 收敛 , 若 不 收敛 , 则 返回 步 又 ai; 

十. 归 一 化 因子 图 像 C, 再 根据 缩放 的 比例 缩放 开 ; 


Of (C,F P 
其 中 ,在 C-step 中 , 令 HCP -0 得 到 
ik 
I, - 
C: F. if 
, 23 kf | Ta 
Ca = ; 
YF, 
f 
在 F-step 中 , 令 ICP =0, 得 到 
LA 
I a 
Y C, : 
a 2 fes) 
Fra IS 
DA 


1.3 唯一 性 约束 

容易 发 现 , 式 1 描述 的 因子 模型 的 结果 在 数学 上 并 
不 是 唯一 的 2 ,因子 和 因子 图 像 的 非 唯一 性 可 以 通过 下 
式 说 明 

I-(CR(R'F)-C'F' (7) 

RrB,CLE oS ECSCIS SD E ARMES EA A] TAB 
阵 ,通过 一 个 可 首 的 和 矩阵 R, 理 论 上 能 得 到 无 数组 符合 
模型 描述 的 因子 图 像 C 和 不 。 因 此 ,我 们 必须 找到 合理 
的 约束 条 件 确保 结果 的 唯一 性 。 

A 传统 方法 :最 小 化 因子 图 像 之 间 的 重合 程度 约束 

sitek.A 假 设 因 子 图 像 之 间 的 重 闭 程度 最 小 (MSO,， 
minimal structure overlap) 的 作为 唯一 性 约束 条 件 ， 
子 间 重 准 程 度 的 计算 公式 为 

e a, |CR| — [CR 


p-lgq-ptli-l Y (cR) Y (CR 
il vw rem i4 


(8) 


FET ACCR, R FREK F RK -T- ES 01 [ROS dE 
负 , 可 用 下 式 表示 : 
f(R)= 2,24 ((CR),) + D 


1 k=1 


0, x20 9) 
ATRUS e x«0 

IgA] ME HER ER Ge bif) , 求 得 R 的 最 优 解 ,并 获 
得 最 终 的 因子 图 像 (CR) 和 代表 每 个 组 织 活 度 曲 线 (RF) 
的 具有 实际 生理 意义 的 因子 。 

B. 新 方法 :基于 动力 学 聚 类 的 约束 

由 于 受到 PET 图 像 分 辨 率 和 部 分 容积 效应 的 影 
响 , 心 肌 信和 号 经 常会 掺 杂 着 10%~15% 的 血液 信号 ,最 小 
化 因子 图 像 之 间 的 重合 程度 约束 未 能 考虑 到 这 一 事 
实 。 因 此 ,本 文 提出 一 种 基于 动力 学 的 时 间 活 度 曲 线 聚 
类 (KB, Kinetic Based) 的 先 验 约束 来 解决 因子 分 析 模 
型 中 的 非 唯一 问题 。 

在 动态 PET 研 究 中 ,每 一 个 像素 的 时 间 行 为 可 以 
通过 TAC 曲线 来 描述 ,因此 ,我 们 基于 像素 的 动力 学 对 
图 像 进行 分 类 。 我 们 的 目标 是 通过 动态 序列 图 像 的 
TAC 曲线 的 形状 和 级 数 对 图 像 进 行 分 类 ,使 得 像素 
TAC 在 同一 类 是 同 质 的 ,而 不 同类 之 间 是 异 质 的 。 假 设 
动态 图 像 1 具 有 KK 个 不 同 特征 的 类 ,模糊 C 均 值 聚 类 算 
法 基于 最 小 化 目标 函数 


| (10) 

其 中 | 代表 欧 氏 距离 ,表示 图 像 7 中 第 j 个 像素 
的 TAC 曲线 ,vy 是 第 k 个 聚 类 的 中 心 的 TAC,， 由 表示 第 7 
个 像素 属于 第 类 的 隶属 度 ,qg(g>1) 是 模糊 参数 。 

目标 函数 的 解 可 以 通过 如 下 迭代 式 得 到 : 

i 设 置 初始 的 类 数 玉 和 停止 准则 s; 

ii. 随 机 初始 化 模糊 隶属 度 和 矩阵 UV 包含 每 一 个 单独 
uj ffs 

iii. 设 置 循环 初始 值 2=-0; 

iv.H U^ T TEXAS UU vi^; 

bU 


Oo. 7 


vR BERE U” 
x 2 N MG - 1) 

eo] =a 

i S 

vi MIR max{0W-UW"}<e 则 停止 ;否则 , 设 b=b+ 
1 返回 第 iv, 

将 聚 类 中 心 w% 通 过 式 (11) 作 为 像素 动力 学 先 验 信 
息 应 用 到 因子 分 析 唯一 性 约束 当中 ， 


fac) DRF) -vill (41) 
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过 最 小 化 目标 函数 (fiict+5j,) 来 解决 因子 分 析 中 
的 非 唯一 问题 , 求 得 R 的 最 优 解 , 最 终 获 得 因子 图 像 
(CR) 和 代表 每 个 组 织 活 度 曲 线 (R"'P) 的 具有 实际 生理 
意义 的 因子 。 


i 


2 实验 

(1) 计 算 机 仿真 数据 

在 计算 机 仿真 实验 中 ,本 文 使 用 基于 如 图 2 所 示 的 
XCAT 数 字体 膜 ,心脏 区 域 由 血 池 ,心肌 组 织 两 个 部 分 
构成 ,在 这 里 将 心脏 区 域 单独 分 离 出 来 。 模 型 的 体 素 点 个 
数 为 64x64x30, 体 素 点 尺寸 大 小 为 3.27 mmx3.27 mmx 
3.27 mm。 我 们 基于 单 房 室 模型 天 模拟 心脏 区 域 生理 代 


AS dH T 


lie 


VE E E sU fs HE RACE HE SR AN "Rb , FE 
241.27 min。 仿真 过 程 中 , 血 输入 函数 以 及 各 组 织 对 应 
的 动力 学 参数 的 测量 值 均 来 自 美国 约翰 霍 普 金 斯 大 学 
PET 医 学 中 心 的 5 位 志愿 者 。 基 于 单 房 室 模型 血 输 
入 函数 以 及 相应 的 动力 学 参数 (设置 动力 学 参数 K= 
1.4822 mL/min/g, k;-0.3159/min 及 血 体积 分 数 为 V; 
0.3829) 生 成 心肌 组 织 的 时 间 活 度 曲 线 , 各 组 织 时 间 活 
度 曲线 及 空间 分 布 如 图 3 所 示 。 整 个 动态 PET 扫 撒 协 
议 设置 为 :12x5 s,6x10 s 共 计 18 个 时 间 帧 。 将 每 一 帧 
血 池 与 心肌 组 织 的 时 间 活 度 赋值 到 其 对 应 组 织 的 空间 
分 布 位 置 上 ,生成 模拟 心脏 区 域 的 动态 PET 图 像 。 

将 动态 图 像 进行 正 向 投影 得 到 随时 间 变 化 的 无 噪 


图 2 (44D XCAT 数 字体 膜 的 正 位 像 ( 左 :男性 , 右 : 女 性 ) 及 (B) 女 性 心脏 区 域 水 平 切片 
Fig.2 Anterior views of the 4D XCAT digital phantom (A; left: male, right: female) and the female 


horizontal plane of myocardial region (B). 


x10 
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图 3 血 池 与 心肌 组 织 的 时 间 活 度 曲线 及 空间 分 布 位 置 ( 横 截面 第 14 层 ) 
Fig.3 Time-activity curve and spatial distribution position of blood pool and 
myocardial tissue (14th slice of the horizontal plane). 


投影 数据 。 投 影 数 据 仿真 过 程 中 产生 的 总 光子 计数 设 
置 为 5x10'。 在 投影 数据 中 加 入 泊 松 噪声 便 可 得 到 仿真 
的 加 噪 投影 数据 。 本 文采 用 FBP 重 建 算法 对 加 品 的 投影 


数据 进行 重建 。 得 到 重建 图 像 之 后 ,用 FWHM=4 mm 
的 高 斯 滤波 器 对 重建 出 的 动态 图 像 进 行 去 噪 。 本 研究 
仿真 实验 均 基 于 该 仿真 数据 。 
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(2) 对 照 方法 及 评价 指标 

在 模型 的 求解 中 的 第 一 步 中 ,引入 了 c 散 度 作 为 测 
度量 化 模型 与 真 值 图 像 之 间 的 距离 ,我 们 将 其 与 传统 的 
最 小 二 乘 测 度 进 行 比 较 。 本 文采 用 归 一 化 误差 2 评价 
不 同 测度 在 每 一 帧 的 剩余 信号 s, 用 归 一 化 误差 2 评价 不 
同 测度 评价 所 有 帧 的 剩余 信号 s。 


(12) 


Nt 

为 了 验证 本 文中 动力 学 聚 类 在 解决 非 唯一 性 问题 
的 准确 性 ,我 们 将 基于 动力 学 聚 类 约束 的 因子 分 析 模型 
与 传统 的 最 小 化 因子 间 重 关 程 度 (MSO) 的 模型 进行 比 
较 。 除 了 视觉 评价 外 ,采用 量化 指标 RMSE 对 因子 与 仿 
真 的 组 织 活 度 曲线 分 布 进行 比较 : 


14.0% [ 
12.094 
..10.096| - 

uU 
8.0% 


6.0% 


4.0% 


Frame 


图 4 不 同 测度 下 模型 在 每 帧 的 归 一 化 误差 % 
Fig.4 Normalization error ð of the model in each 
frame under different measures. 


3.2 唯一 性 约束 

得 到 初始 因子 后 ,用 最 小 化 因子 间 重 县 程度 
(MSO) 和 基于 动力 学 聚 类 (KB) 的 两 种 约束 方法 解决 因 
子 分 析 模 型 中 的 非 唯 一 问题 。 图 6、 图 7 分 别 给 出 了 对 
我 也 燥 比 (SNRJ)=10 的 动态 序列 用 两 种 不 同 约束 方法 
(模型 的 测度 为 a 散 度 测度 , 且 a=25) 所 得 到 的 因子 和 
子 图 像 ,由 图 可 知 ,所 得 出 的 因子 1 与 因子 2 分 别 对 应 于 
血 池 TAC 与 心肌 TAC, 所 得 到 的 因子 图 像 1 与 因子 图 像 
2( 图 中 绘制 的 是 第 14 层 横断 面 图 ) 也 对 应 于 心脏 的 血 
和 心肌 组 织 , 对 于 血 池 的 活 度 曲 线 , 两 种 方法 所 得 到 
因子 的 精确 程度 无 显著 差别 ,而 对 于 心肌 活 度 曲 
线 , 用 KB 约束 得 到 的 因子 比 MSO 约束 得 到 的 因子 精 
确 度 高 。 


IL. 2 
cmo 


y] norm(F)-norm(F, e 
RMSE = < 
Y Wnorm(F ur) 
其 中 F 代 表 因子 曲线 ,Ps 代表 各 组 织 的 活 度 曲线 。 
WV g— WRR TRE RT , norm CF NAE F ETTIH— 
化 , 即 norm(F) SF . 


x100% (13) 


J 


3 结果 
3.1 o 散 度 与 最 小 二 乘 测度 建 模 的 精确 性 的 比较 

通过 最 小 化 模型 与 真实 图 像 之 间 的 测度 得 到 初始 
因子 和 因子 图 像 后 ,可 用 剩余 信和 号 来 评价 不 同 测度 的 优 
劣 。 图 4 是 最 小 二 乘 测度 与 a 散 度 测度 (a=1,10,20,30， 
40,50) 所 得 出 的 模型 归 一 化 误差 6,, 图 5 是 动态 图 像 所 
有 帧 的 归 一 化 误差 6, 随 oe 的 变化 情况 ,由 图 可 知 , 仪 当 o= 
1 时 a 测度 与 最 小 二 乘 测度 的 归 一 化 误差 相等 , 当 a 取 其 
它 值 时 ,a 散 度 测度 的 结果 的 归 一 化 误差 均 优 于 最 小 二 
乘 测度 。 


8.50% f 
— a divergence measure 


8.40% f — Least square measure 
8.30% 上 
8.20% 上 
8.1095 上 
8.00% 上 


7.9096 上 


7.8096 
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a 


图 5 不 同 测度 下 模型 的 归 一 化 误差 8 
Fig. Normalization error à; of the model under 
different measures. 


3.3 a 值 的 选取 

为 了 进一步 验证 a 散 度 在 本 模型 中 的 有 效 性 基 选 取 
最 优 的 a 值 ,本 文选 取 不 同 的 a 值 在 两 种 约束 模型 下 进行 
因子 分 析 , 为 避免 实验 受 噪声 的 随机 性 影响 ,针对 每 种 
方法 运行 30 次 求 平均 值 ,绘制 出 不 同 w 值 在 MSO 和 KB 
两 种 约束 下 各 因子 的 RMSE 图 8, 在 KB 约束 中 ,心肌 因 
子 的 RMSE 基本 保持 稳定 , 当 o=5 时 估计 的 血 池 因子 
RMSE 最 小 ,因此 选取 5 作为 本 文 的 最 优 的 a 值 ;而 在 
MSO 约 束 中 ,两 因子 在 o>12 时 基本 趋 于 稳定 。 选 取 了 
最 优 的 值 后 ,将 a 散 度 测度 与 最 小 二 乘 测度 在 两 种 不 同 
约束 下 作 比 较 见 表 1, 结 果 说 明 a 散 度 的 RMSE 比 最 小 
二 乘 的 低 , 即 该 测度 更 能 提取 出 更 准确 的 结果 。 此 外 ， 
由 图 8 可 知 无 论 c 取 何 值 ,用 KB 约束 模型 得 出 的 两 因子 
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图 6 用 MSO 约束 所 得 因子 及 因子 图 像 
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Fig.6 The factors and factor images obtained by MSO constraints. 
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图 7 用 KB 约束 所 得 因子 及 因子 图 像 


Fig.7 Factors and factor images obtained by KB constraints. 


的 平均 RMSE 曲线 位 于 MSO 的 下 方 , 这 亦 说 明了 动力 
学 聚 类 约束 要 优 于 小 化 因子 间 重 闭 程 度 约束 。 
3.4 模型 对 噪声 的 敏感 度 

为 了 评价 噪声 对 模型 的 影响 ,我 们 通过 在 图 像 的 投 
影 域 中 加 入 不 同 规模 的 噪声 ,重建 出 不 同 信 噪 比 的 图 
像 。 表 2 表示 对 不 同 信 噪 比 的 动态 序列 ,用 MSO 和 KB 
两 种 不 同 约束 所 得 到 的 各 因子 的 RMSE。 由 表 可 知 , 当 
噪声 逐渐 增 大 时 ,因子 的 RMSE 也 相应 增 大 ,估计 的 血 
池 的 时 间 活 度 曲 线 比 心肌 活 度 曲线 精确 。 


4 讨论 

a 散 度 测度 的 目的 是 将 动态 图 像 做 一 个 初步 的 分 
解 ,由 于 分 解 的 结果 具有 非 唯 一 性 ,所 以 根据 聚 类 中 心 
作为 先 验 信息 将 其 中 一 组 解 线 性 变换 为 具有 实际 生理 


不 能 作为 最 终 的 TAC ,因为 聚 类 与 聚 类 中 心 过 于 粗糙 ， 
不 能 较 精 确 地 反应 真实 的 动态 图 像 ,尤其 是 对 于 噪声 较 
大 的 动态 图 像 。 

在 加 入 唯一 性 约束 时 ,我 们 引入 一 个 惩罚 项 来 保 
证 我 们 得 到 的 结果 是 非 负 的 ,对 于 两 个 正则 化 参数 b, 和 
bs, 实验 发 现 当 b,,bye[10, 500] 时 ,结果 无 明显 差异 ,鉴于 
此 在 本 实验 中 将 b,,b; 均 设置 为 100。 对 于 不 同 的 数据 ， 
我 们 可 以 参考 本 数据 选取 合适 的 bp 值 ,大 结果 出 现 负 值 ， 
我 们 只 需要 适当 地 增加 2 的 值 保证 结果 是 非 负 的 即 可 。 

从 该 心脏 体 膜 中 提取 的 结果 我 们 我 发 现 ,在 同样 规 
模 噪 声 的 情况 下 ,估计 出 的 血 池 的 时 间 活 度 曲线 比 心肌 
组 织 的 活 度 曲线 要 精确 ,主要 原因 是 血 池 的 体 素 数 比 心 
肌体 素数 多 , 受 部 分 容积 效应 及 噪声 的 影响 较 小 。 在 心 
肌 壁 的 某 些 部 分 仅 有 2~3 个 体 素 的 厚度 , 极 容易 受到 部 


意义 的 时 间 活 度 曲 线 。 值 得 注意 的 是 ,动力 学 聚 类 绪 


分 容积 效应 及 噪声 的 影响 降低 结果 的 精确 度 。 
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图 8 不 同 o 值 在 MSO 和 KB 两 种 约束 下 各 因子 的 RMSE 
Fig.8 RMSE of the obtained factors with MSO constraints and KB 
constraints undera divergence measure. 
表 1 两 种 测度 下 应 用 MSO 约束 和 KB 约束 的 所 得 各 因子 的 。 准确 性 有 积极 的 作用 ,在 高 噪声 下 (SNR=5) 仍 能 准确 的 


RMSE 
Tab.1 RMSE of the obtained factors with MSO constraints and KB 
constraints under two measures 


Least square measure adivergence measure (0=5) 


Item 

MSO KB MSO KB 
Blood 4.0796 1.696 2.4396 0.52% 
Myocardium 4.89% 1.48% 4.43% 1.45% 


表 2 不 同 信 噪 比 下 各 因子 的 RMSE 
Tab.2 RMSE of various factors under different SNRs 


No noise SNR-10 SNR-5 
Item 
MSO KB MSO KB MSO KB 
Blood 0 0 2.9090 0.64% 8.96% 1.3396 
Myocardium 0 0 44499 1.44% 11.32% 0.467946 


Notes: MSO representsMinimal Spatial Overlap constraint method and 
KB represents Kinetic Basedconstraint method. 


本 研究 提出 一 种 基于 动力 学 聚 类 与 o 散 度 测度 的 动 
三 维 心 脏 PET 图 像 因 子 分 析 研 究 模 型 ,实现 了 对 心脏 
PET 图 像 各 组 织 的 时 间 活 度 曲 线 进 行 精确 提取 。 传 统 
的 因子 分 析 模 型 利用 最 小 二 乘 模型 ,用 最 小 化 因子 图 像 
dvi uus ue UN NM ,本 

过 最 小 化 动态 图 像 与 因子 分 析 模 型 的 wx 散 度 测度 ， 
ae 同 的 c 值 选取 比 最 小 二 乘 适 合 的 测度 得 到 初 
始 因子 及 因子 图 像 ,并 在 旋转 因子 时 引入 PET 图 像 像素 
的 动力 学 聚 类 先 验 信 息 ,我 们 将 模型 应 用 于 基于 XCAT 
的 仿真 数据 上 ,通过 实验 分 析 发 现 ,用 a 散 度 测度 以 及 基 
于 动力 学 聚 类 的 先 验 约束 均 对 组 织 活 度 曲 线 的 提取 的 


态 


提取 出 各 组 织 的 时 间 活 度 曲线 ,用 视觉 评价 及 量化 评价 
等 方法 优势 明显 。 然 而 ,本 方法 在 旋转 因子 时 引入 正则 
化 参数 无 法 自 适应 选择 ,在 下 一 步 工 作 中 将 进一步 研究 


正则 化 参数 自 适 应 选择 问题 。 
参考 文献 : 
[1] 王 荣 福 PET/CT: 分 子 影像 学 新 技术 应 用 [M] 北京: 北京 大 学 医学 出 


版 社 , 2011: 108-13. 
[2 


Lu 


Su Y, Arbelaez AM, Benzinger TLS, et al. Noninvasive estimation 
of the arterial input function in positron emission tomography 
imaging of cerebral blood flow [J]. J Cereb Blood Flow Metab, 
2013, 33(1): 115-21. 

[3] Schupbach WMM, Emese B, Loretan P, et al. Non-invasive 
diagnosis of coronary artery disease using cardiogoniometry 
performed at rest|J]. Swiss Med. Wkly, 2008, 138(15-16): 230-8. 

[4] Bisker J. Principles and practice of positron emission tomography 
[J].Am J Roentgenol, 2012, 180(5): 1238. 


[5] Garbarino S, Caviglia G, Brignone M, et al. Compartmental analysis of 
renal physiology using nuclear medicine data and statistical 
optimization|J ]. Physics, 2012, 3(24): 625-37. 

[6] Fakhri GE, Sitek A, Guerin B, et al. Quantitative dynamic cardiac 
82Rb PET using generalized factor and compartment analyses[J]. J 
Nucl Med, 2005, 46(8): 1264-71. 

[7] Moody J, Lee BC, Corbett JR, et al. Precision and accuracy of 


clinical quantification of myocardial blood flow by dynamic PET: A 
technical perspective[J]. J Nucl Cardiol, 2015, 22(5): 935-51. 
[8] Tio RA, Dabeshlim A, Siebelink HM, et al. Comparison between 


the prognostic value of left ventricular function and myocardial 


perfusion reserve in patients with ischemic heart disease|J]. J Nucl 
Med, 2009, 50(2): 214-9. 

[9] Dehak N, Kenny PJ, Dehak R, et al. Front-end factor analysis for 
speaker verification[J]. IEEE Trans Speech Audi P, 2011, 19(4):788- 
98. 

[10] Kim J, Herrero P, Sharp T, et al. Minimally invasive method of 


"i 


201801.00603v1 


chinaXiv 


: 1584 - 


J South Med Univ, 2017, 37(12): 1577-1584 


ChinaXiv 合 作 期 刊 


http://www.j-smu.com 


determining blood input function from PET images in rodents[J]. J 
Nucl Med, 2006, 47(2): 330-6. 

[11] Klein R, Beanlands RS, Wassenaar RW, et al. Kinetic model-based 
factor analysis of dynamic sequences for 82-rubidium cardiac 
positron emission tomography [J). Med Phys, 2010, 37(8): 3995- 
4010. 

[12] Ketata I, Sallemi L, Morainnicolier F, et al. Factor analysis-based 
approach for early uptake automatic quantification of breast cancer 
by I8F-FDG PET images sequence [J]. Biomed Signal Proces, 
2014: 19-31. 

[13] Paola RD, Bazin JP, Aubry F, et al. Handling of dynamic sequences 

in nuclear medicine| J]. IEEE Trans Nucl Sci, 1982, 29(4): 1310-21. 

[14] Wu HM, Hoh CK, Choi Y, et al. Factor analysis for extraction of 

blood time-activity curves in dynamic FDG-PET studies[J]. J Nucl 

Med, 1995, 36(9): 1714-22. 

[15] Whittle P. On principal components and least square methods of 

factor analysis|J]. Scand Actuar J, 2012, 35(3-4): 223-39. 


[16] Lord S, Galna B, Verghese J, et al. Independent domains of gait in 
older adults and associated motor and nonmotor attributes: 
validation of a factor analysis approach|J). J Gerontol Ser A, 2013, 
68(7): 820-7. 

[17] Mabrouk R, Dubeau F, Bentabet L. Dynamic cardiac PET imaging: 
extraction of time-activity curves using ICA and a generalized 
Gaussian distribution model[J]. IEEE Trans Biomed Eng. 2013, 60 
(1): 63-71. 

[18] Su Y, Welch M J, Shoghi K, et al. The application of maximum 

likelihood factor analysis (MLFA) with uniqueness constraints on 

dynamic cardiac microPET data [J]. Phys Med Biol, 2007, 52(8): 

2313-34. 


[19] Sitek A, Gullberg GT, Huesman RH. Correction for ambiguous 
solutions in factor analysis using a penalized least squares objective. 
[J]. IEEE Trans Med Imaging, 2002, 21(3): 216-25. 

[20] Boutchko R, Mitra D, Baker SL, et al. Clustering-initiated factor 
analysis application for tissue classification in dynamic brain 
positron emission tomography [J]. J Cereb Blood Flow Metab, 
2015, 35(7): 1104-11. 

[21] Jia S, Qian Y. Constrained nonnegative matrix factorization for 
hyperspectral unmixing|J|. IEEE Trans Geosci Electron, 2009, 47 
(1): 161-73. 

[22] Critchley F, Marriott P. Computing with Fisher geodesics and 
extended exponential families[J]. Stat Comput, 2014, 26(1): 1-8. 
[23] Teng Y, Zhang T. Generalized EM-type reconstruction algorithms 
for emission tomography [J). IEEE Trans Med Imaging, 2012, 31 

(9): 1724-33. 

[24] Lu L, Karakatsanis NA, Tang J, et al. 3.5D dynamic PET image 

reconstruction incorporating kinetics-based clusters [J]. Phys Med 

Biol, 2012, 57(15): 5035-55. 

[25] Segars WP, Mahesh M, Beck TJ, et al. Realistic CT simulation 

using the 4D XCAT phantom[ J]. Med Phys, 2008, 35(8): 3800-8. 

] 马 晓 勉 , 路 利 军 , 高 园 园 ,等 . 基于 低 秩 框架 的 心肌 灌注 动态 PET 图 像 

恢复 [可 . 暨南 大 学 学 报 : 自然 科学 与 医学 版 , 2016, 37(1): 78-86. 

[27] El Fakhri G, Trott CM, Sitek A, et al. Dual-tracer PET using 


[2 


CN 


generalized factor analysis of dynamic sequences |J]. Mol Imaging 
Biol, 2013, 15(6): 666-74. 

[28] Watabe H, Ikoma Y, Kimura Y, et al. PET kinetic analysis-- 
compartmental model[J]. Ann Nucl Med, 2006, 20(9): 583-8. 


(编辑 : 吴 锦 雅 ) 


